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1. Introduction 

The problem of quantization of gravitational field is one of the fundamental problems 
in theoretical physics. One of the approaches to the construction of a theory of 
quantum gravity is called Loop Quantum Gravity (LQG), which began with the 
seminal work by Ashtekar [1] and has since developed in a large research discipline 
[2]. Research in LQG is focused on two complementary fronts. The approach used 
on one front is to perform the canonical quantization of the gravitational field, which 
is known as canonical LQG. This is done by rewriting General Relativity (GR) in 
the Hamiltonian form, by using connection variables and their momenta, with the 
goal of constructing the appropriate physical Hilbert space. This amounts to working 
with the the spin-network states. On the other front, the approach is to perform the 
path-integral quantization of GR by using the connection variables, and this requires 
using a spacetime triangulation where the triangles carry spin labels [3] . The resulting 
theories are called spin-foam models. 

The spin-foam path integrals can be viewed as transition amplitudes from one 
spin-network state to another one, via the corresponding spin foam. This picture 
provides a relationship between the canonical and the path-integral formulations. 
Although one can couple fermions in canonical LQG, it is impossible to couple correctly 
the fermions in the spin-foam formulation. This is a consequence of the fact that there 
are no edge lengths in a spin-foam model, which is related to the absence of the tetrads 
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in the Plebanski action. Furthermore, the classical limit of a spin-foam model is the 
area-Regge theory HIS], which is based on a twisted geometry, where metric is not 
always defined. 

Recently a new idea has been proposed for the spin-foam approach, in the form 
of the spin-cube models [51E]. These models represent a categorical generalization of 
the spin-foam models. One can generalize the path integral construction based on a 
group to the construction based on a 2-group. For the basic notions of 2-group theory 
see mi- The construction of the spin-cube path integral rests on a reformulation of 
GR as a constrained BFCG theory. The BFCG theory is a categorical generalization 
of the well-known BF theory. The consequence for the path integral is that it becomes 
a sum over the labels for the edges, triangles and the tetrahedrons in a triangulation. 
Hence one obtains a sum of amplitudes for a colored 3-complex, therefore the name 
spin cube. This also implies that a spin-cube model can describe a transition amplitude 
from a spin-foam state to another spin-foam state, so that one would like to study a 
canonical quantization of the BFCG formulation of GR in a spin-foam basis. 

In particular, given a Lie group G and its Lie algebra g, one can assign to it an 
action of the form 

Sbf = f {B A F)g , (1) 

Jm 

where H is a g-valued Lagrange multiplier two-form, while F = dA + H A H is the 
curvature two-form for the g-valued connection one-form A. The ( , )g represents 
the invariant nondegenerate symmetric bilinear form in g. The BF theory relevant 
for the construction of spin-foam models is based on the Lorentz group 50(3,1). A 
categorical generalization of the BF theory is based on the concept of a 2-group, 
constructed of two groups Q and H in a particular way (see [7] for details). It is called 
the BFCG theory [9l[T0], and has the action 

Sbfcg = [ {B A F)g {C AG)t, ■ (2) 

Jm 

Here the second term consists of a l)-valued one-form Lagrange multiplier C, and a 
curvature three-form G — dfl A A ft for the t)-valued two-form /3, where t) is the 
Lie algebra of the group H. The pair {A, /3) is called the 2-connection of the given 
2-group, while the pair {F,G) is the corresponding 2-curvature. The ( , )[, is the 
invariant nondegenerate symmetric bilinear form in f), which is g-invariant. 

The importance of the BFCG theory for the Poincare 2-group, defined by the 
choice G = 50(3,1) and 77 = K"*, lies in the fact that one can construct the action for 
GR by simply adding an additional term to the BFCG action, called the simplicity 
constraint: 

5=/ (HAi?)g + (eAG)f,-((^A(H-*(eAe)))g. (3) 

Jm 

Here we have relabeled C = e and F = R, since in the case of the Poincare 2-group 
these fields have the interpretation of the tetrad field and the curvature two-form for 
the spin connection A = uj. The g-valued two-form (f is an additional Lagrange 
multiplier, featuring in the simplicity constraint term. The * is the Hodge dual 
operator for the Minkowski space. See [6] for details. 

The construction of ([3]) is in full analogy to the Plebanski model, where GR is 
constructed by adding a suitable simplicity constraint to the BF theory based on 
the Lorentz group. However, in contrast to the Plebanski model, the constrained 
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BFCG model has one big advantage. Namely, the Lagrange multiplier C has the 
interpretation of the tetrad field, which is therefore explicitly present in the topological 
sector of the action. This is not the case for the Plebanski model, where the 
simplicity constraint merely infers the implicit existence of the tetrad fields. Upon 
the covariant quantization, the Plebanski action gives rise to spin-foam models, while 
the constrained BFCG action gives rise to the spin-cube model. Then, the explicit 
presence of the tetrad in ([3]) enables us to easily couple matter fields to gravity in 
the spin-cube model, in contrast to spin-foam models where this is a notoriously hard 
problem. 

As a classical theory, the constrained BFCG action lends itself also to the 
canonical quantization programme. In the canonical approach, the first and crucial 
step is to perform the Hamiltonian analysis of the theory, study the algebra of 
constraints, and eliminate second class constraints from the theory. However, due 
to the technical complexity of the Hamiltonian analysis, it is wise to discuss the pure 
BFCG theory first, leaving the constrained theory for later. That is the aim of this 
paper. We will perform the full Hamiltonian analysis of the unconstrained BFCG 
theory based on the Poincare 2-group, as a preparation for the more complicated case 
of the constrained BFCG theory ([3]). The similar analysis has been done for the 
BF theory in and our analysis represents the generalization of that work to the 
BFCG case. 

We should note that the Hamiltonian analysis of the BFCG model has already 
been done in a gauge-fixed form in [12]. In this paper we improve those results by 
providing a gauge-invariant canonical analysis. 

There is an interesting relationship between the BFCG theory for the Poincare 
2-group on one hand, and the topological Poincare gauge theory on the other. Perhaps 
surprisingly, the two theories are equivalent, while their Hamiltonian structure is vastly 
different. This was discussed to an extent in [12] . but the full Hamiltonian analysis 
presented here lends itself nicely to a more complete comparison of the Hamiltonian 
formulations for the BFCG theory and the topological Poincare gauge theory. The 
results of this comparison are especially intriguing, and provide additional insight into 
the structure of the theory. 

The paper is organized as follows. In section [2] we introduce in detail the actions 
for the BFCG theory and the topological Poincare gauge theory, give a short overview 
of the Lagrange equations of motion, and prepare for the Hamiltonian analysis. The 
bulk of the Hamiltonian analysis is done in section [Sj We evaluate the conjugate 
momenta for the fields, obtain the primary constraints and construct the Hamiltonian 
of the theory. Then we impose consistency conditions on the primary constraints, 
which leads to secondary constraints and some determined Lagrange multipliers. The 
consistency conditions of the secondary constraints turn out to be satisfied identically, 
and do not introduce any new constraints. The constraints are then separated into 
first and second class, and their algebra is computed. Finally, the number of physical 
degrees of freedom is calculated, and ends up being zero, confirming that the theory 
is indeed topological. Building on these results, in section |4| we construct the Dirac 
brackets, which facilitate the elimination of the second class constraints from the 
theory, leading to the reduction of the phase space. Section [5] is devoted to the study 
of the properties of the reduced phase space, with the emphasis on the differences 
between the BFCG model and the topological Poincare gauge theory. Section [51 
contains the discussion of the results and our concluding remarks. The Appendix 
provides some technical details about the derivation and the discussion of the Bianchi 
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identities, used in the main text. 

Our notation is as follows. The spacetime indices are denoted with lowercase 
Greek alphabet letters pL,v,p,... and take the values 0,1, 2, 3. When discussing the 
foliation of spacetime into space and time, the spacetime indices are split as /r = (0, i), 
where the lowercase indices from the middle of the Latin alphabet i,j, fc,... take only 
spacelike values 1,2,3. The Poincare group indices are denoted with lowercase letters 
from the beginning of the Latin alphabet, a,b,c,... and take the values 0,1, 2, 3. They 
are raised and lowered with the Minkowski metric rjab = diag(—1,1,1,1). Capital 
Latin indices A, B,C,... represent multi-index notation, and are used to count the 
second class constraints. Antisymmetrization is denoted with the square brackets 
around the indices, 


1 


-^[ab] rj {^ab Aba) ■ 


(4) 


In order to simplify the notation involving Poisson brackets, we will adopt the following 
convention. The left quantity in every Poisson bracket is assumed to be evaluated 
at the point x = {t,x), while the right quantity at the point x' = {t,x'). In 
addition, we use the shorthand notation for the 3-dimensional Dirac delta function 
^(3) = ^(3) (a; _ x'). This allows us to write an explicit but bulky expression like 

{ A{t, x ), B{t, x') } = C(t, {x — x') (5) 

more compactly as 

{A,B} = C5^^\ ( 6 ) 

without ambiguity. This notation will be used systematically unless explicitly stated 
otherwise. 

2. BFCG action 


The BFCG theory for the Poincare 2-group is defined by the action 


Sbfcg= / Bab^R‘^^ + e^ ^Ga. (7) 

JM 

The variables of this action are the one-forms e“, and the two-forms 5“^, /3“. The 
curvatures R°‘^ and G“ are the field strengths of the 2-connection (a;“^,/3“), 

= dcu“'' -b a;“e A , (8) 

= d/ 3 “ -b a;“b A . ( 9 ) 

The fields and e“ play the role of the Lagrange multipliers. Usually one would 
denote the latter multiplier as G“, but we shall instead label it as e“ since it will be 
interpreted as the tetrad field. Similarly, the usual notation for the connection one- 
form and its field strength is A and F respectively, but in our case they are denoted 
and i?“^, since they are interpreted as the spin connection and the Riemann curvature 
two-form. 

It is also convenient to introduce torsion as the field strength for the tetrad e“, 

pa ^ yga ^ ^ ^ ^b 

Then, performing a partial integration in the second term in © and using the Stokes 
theorem one can rewrite the action as 

ab 


StPGT = / Bab A R -b / 3 “ A Ta — 


e“ A/3a 


( 11 ) 
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where now and play the role of Lagrange multipliers. Aside from the immaterial 
boundary term, this action represents the topological Poincare gauge theory (TPGT). 
In order to fully appreciate the relationship between the two theories in the sense of 
the Hamiltonian analysis, let us introduce a parameter ^ G R and rewrite the action 
as 


S= j Sab A + ee“ A Ga + (1 - 0/3“ A Ta , (12) 

Jm 

where we have dropped the boundary term. It is obvious that the action (HU is a 
convenient interpolation between 0 and (ED, to which it reduces for the choices ^ = 1 
and ^ = 0, respectively. The action ED will therefore be the starting point for the 
Hamiltonian analysis. 

It is also clear that all three actions o, ED and ED give rise to the same set of 
equations of motion, since these do not depend on the boundary. Taking the variation 
of (TT^ with respect to all the variables, one obtains 


where 


5B ■. = 0, 

(5/3 : r“ = 0 , 

5e: G“ = 0, 

5uj : - e[“ A /3'’1 = 0 , 


= dB“^ + A A . 


(13) 


(14) 


The first two equations of motion tell us that spacetime has Minkowski geometry, 
while the third and fourth determine /3“ and B°‘^. As we shall see later, there are no 
local propagating degrees of freedom in the theory. 

Finally, for the convenience of the Hamiltonian analysis, we need to rewrite both 
the action and the equations of motion in a local coordinate frame. Choosing as 
basis one-forms, we can expand the fields in the standard fashion: 

e“ = = w^^dx'^, (15) 

^ A dx" , /3“ = i/SV^dx'^ A dx" . (16) 


Similarly, the field strengths for w, e and /3 are 

^dx'', 

T“ = ^TVdx'^Adx", 

G“ = iG“^;ypdx^ A dx"" A dx^ . 

Using the relations ([8]), ([9]) and ED, we can write the component equations 

TDah _ , ,ab , ,ab i , ,a , ,cb , ,a , ,cb 

^ — O^UJ y — OyUJ ^ ly — id ^ 1 

-b - w“b,ye'’p , 

G“pi/p = dfj_f5°‘yp -b di,/3°‘pp -b dpf5°‘pi, -b vp + bo^^bvfi^pp + ^°‘bpl3^pu ■ 

Substituting expansions ED, (ED and ED into the action, we obtain 


5 = 




IM 


-^BabpvR°'^ p(T + ^S.apG°‘ i,prj H - 4”^/^' 


T a 

p<T 


(17) 


(18) 


(19) 
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Assuming that the spacetime manifold has the topology A4 = S x R, where E is 
a 3-dimensional spacelike hypersurface, from the above action we can read off the 
Lagrangian, which is the integral of the Lagrangian density over the hypersurface E: 


L = 





ab^u 


j^ab 


pa 


-p (7“ 

vpa 


1-^ 

4 



( 20 ) 


Finally, the component form of equations of motion is: 


= 0, 


rpa _ n 
^ pu - 5 


— 0 

pup - 5 


^\pup 




= 0 . 


( 21 ) 


3. Hamiltonian analysis 


Now we turn to the Hamiltonian analysis of the BFCG theory. A detailed review of 
the general formalism can be found for example in |13| . Chapter V. In addition, the 
equivalent procedure for the ordinary BF theory has been done in m- 

As a first step, we calculate the momenta tt corresponding to the field variables 
e“^, 0 ;“^^ and Differentiating the Lagrangian with respect to the time 

derivative of the appropriate fields, we obtain the momenta as follows: 


7r{B)abP'' 

'^{e)aP 

T^{^)abP‘ 


SL 

6L 

6doe°-)_c 

5L 

SL 


0 , 

2 ^ Havp 1 

e°P''PBab.p , 
-ie^P'^Peap ■ 


( 22 ) 


None of the momenta can be solved for the corresponding “velocities”, so they all give 
rise to primary constraints: 


P{B)abP‘' 


P{e)aP 

= <e)aP-^^e°P''PPa.p-d 

P{co)abP 

= TT{u:)abP - e^P’^PBabup ~ 0 , 

P^aP'' 

= AP)aP'' +ie°P’'Peap^0. 


The weak, on-shell equality is denoted “ss”, as opposed to the strong, off-shell equality 
which is denoted by the usual symbol “=”. 

Next we introduce the fundamental simultaneous Poisson brackets between the 
fields and their conjugate momenta. 


{ . <B)cdP^ } 


, 7r(e)6'^ } 


Tr{oj)cd''} 





( 24 ) 
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and we employ them to calculate the algebra of primary constraints, 

{ P(e)“'=, P(/3)&b } = 

while all other Poisson brackets vanish. Note that the algebra of primary constraints 
is independent of 


(25) 


Next we construct the canonical, on-shell Hamiltonian: 

'1 


Hr^ = 


d^a; 


3- 7r(i?),b'^"aoi?“''^,+7r(e)a'^5oe%+ 

JY. 


(26) 


-L. 


The factors 1 /4 and 1/2 are introduced to prevent overcounting of variables. Using (IT^ 
and (|20|. one can rearrange the expressions such that all velocities are multiplied by 
primary constraints, and therefore vanish from the Hamiltonian. After some algebra, 
the resulting expression can be written as 


Hr = - d3a;e°bfc 


-^BabOiR^^^jk + geaoG°ijfe + 


+ —l3a0kT°'ij + —l^abO {'^iB°‘^jk — e°'if3^jk) 


(27) 


up to a boundary term. The canonical Hamiltonian does not depend on any momenta, 
but only on fields and their spatial derivatives. Also, note that it does not depend on ^ 
either. Finally, introducing Lagrange multipliers A for each of the primary constraints, 
we construct the total, off-shell Hamiltonian: 

1 

2 ' 

1 


Ht = Hc-\- J d^x 

1 

1 


A(e)%P(e)a^ + -A(ujr\P(u;)ab^ 


+ ^A(H)“V,P(B)ab^" + lA(/3)%,P(/3)a'^" 


(28) 


We proceed with the calculation of the consistency requirements for the primary 
constraints, 

P = {P, Hrj-O. (29) 

Half of the consistency requirements will give the secondary constraints S, while the 
other half will determine some of the multipliers A. In particular, requiring that 

P(e)a°^0, 

)ab^ 

we obtain the following secondary constraints: 


PiB)ab°^^0, 

P(/3)„°*«0, P{uj)ab°-0, 


(30) 


5(P)“V = P“V«iO, 


0 , 


= ^ 0 , 


SiGY 

SiT)\j 

S{BY'^ = eObfc « 0. 

The remaining consistency conditions for the primary constraints, 

PiB)ab^'^ - 0 , P(e)a'= 0 , 

P(/3)a^'= 0 , P(cc),/ « 0 , 


(31) 


( 32 ) 
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\{e)\ 


2V[,rz]o - , 

Vie“o - w“boe^ , 

2V[jS“^]o + ■ 


This leaves the multipliers 

A(/3)“o,, A(ero, X{BT\^, (34) 


undetermined. 

As the next step we impose the consistency conditions for the secondary 
constraints dSH), 


5(r)“y«o, s{BY^KiQ. 


After a straightforward but lengthy calculation, it turns out that all these conditions 
are identically satisfied, producing no new constraints and determining no additional 
multipliers. Therefore, at this point all the consistency conditions have been 
exhausted. 

Once we have found all the constraints in the theory, we need to classify them 
into first and second class. While some of the second class constraints can already be 
read from (ESD, the classification is not easy since constraints are unique only up to 
linear combinations. The most efficient way to tabulate all first class constraints is to 
substitute all determined multipliers into the total Hamiltonian (1281) and rewrite it in 
the form 


Hf 



\\{BY\^ c^{B)ab^ + A(e)“o f>{e)a + A(/3)“i cf[l3)Y 
+ \\{ioY^ - \BaW^ - eao 

-l3a0^f^iTr -^UJawH^BY^ . 


(36) 


The quantities (f are linear combinations of constraints, but must all be first class, 
since the total Hamiltonian weakly commutes with all constraints. Written in terms 
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= P{B)alP^ , 

4>{e)a = P{e)a ^, 

= P{P)a^^ , 

(j){uj)ab = P{i^)ab° , 

HGT = + V,P(e)“ - , 

6 4 

HTr = i£°*^'=^(r)V-V,P(/3)“^ + ^e6,P(P)“^*^ 

(j){S/BY^ = 5'(P)“*' +V*P(u;)“'’* 

-2e[“,P(e)'’l* -/3[“yP(/3)*'l*^'. 

These are the first class constraints in the theory. The remaining constraints are 
second class: 

x{B)ab^^ = P{B)ab ^'^, x(e)a* = P{e)G , 

(38) 

X{oj)aG = P{co)ab^ , X(/3)a*^' = P(/3)a*^' ■ 

In order to calculate the full algebra of constraints, it is convenient to express them 
as functions of fundamental variables, as follows: 


YB)aG = 7r{B)ab°^ , 

(l){e)a = 7r(e)a°, 

<('(/3)a* = 7r(/3)a°* , 

'('(w)a& = 7I'(w)a6° , 

(j){RY>>i = - Vj7r(P)“^*^' 


4>{GY 

(t){TY" 


= ^Y-i^GYjk + V,7r(e)“ - i/36,,7r(P)“*'*^' , 

= - V,7r(/3)“^' + ie6,7r(P)“'’*^ 


(j){\7BY^ = V,7r(a;)“^* - p[“cy7r(P)'''“^' - 2e[“i7r(e)'’l* - /3[“i,-7r(/3)''l*^' 


X{B)ab^’^ 

= HbyyY 

x(e)Y 

= ^(e)a* - 

x(^)ab" 

= 7r{u;)aY - , 

X(/3)a"^ 

= ^(/3)a*^+Ce°*^'"eafe. 
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The algebra between the first class constraints is then 


{cf{TY\ (fiVBU) 
{4>{RY'>\ct>{VBU} 

{ cf{VBY^ , } 


-45[Xi?)''l,f5(3) , 
-45[XVB)''ld]5(3), 


the algebra between the second class constraints is, according to (|25]) . 


(41) 






between the first and second class constraints is 

{cfiRr»\xibo)cF} 

= 45[:x(i3)''ld]*'d("(, 

{cf{Gr,xiboU} 

= 25“x(e)d]*5(3), 

{<^(G)“,x(/3)c^"} 

= -ix(i?) V"d(3), 

{cf{Tr\x{oo)cd^} 

= -2S^^xWdY6^^\ 

{cfiTYGxieW} 

= ix(S)“b*^d(3), 

{(^(Vi?)“^ x(w)c/} 

= 45[);x(ca)rf]''l*5(3) , 

{ cfiVBY^ , x(/3)c^" } 

= -25i“x(^)'''^'''d(3) , 

{<^(Vi?)“^x(e)e*} 

= -25i“x(e)'''*5(3), 

{ ^{VBY^ , x(S)cd^" } 

= 45[:x(S)d]''''"d(3). 


(42) 


(43) 


All other Poisson brackets among (f and x zero. 

We see that the algebra is closed, and all Poisson brackets involving (f constraints 
weakly vanish, confirming that all (j) are indeed first class. Also, the Poisson brackets 
between y constraints do not weakly vanish, confirming that x ^^e indeed second class. 
Finally, note that the structure constants do not depend on Y despite the fact that 
the constraints (f and x do. 

The last main step in the Hamiltonian analysis is the counting of the physical 
degrees of freedom. Given N initial independent fields in the theory, the dimension of 
the full phase space is 2N. From this one subtracts the total number F of first class 
constraints, the total number S of second class constraints, and the total number F 
of gauge fixing conditions. The result is the dimension of the physical phase space, 
2n, where n is the number of physical degrees of freedom. Thus we have the general 
formula (see for example m). 

n = N - F - - . 

2 


(44) 
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The number of independent field components for each of the fundamental fields is 




pd 

Tjab 

24 

24 

16 

36 


which gives the total N = 100. Similarly, the number of independent components for 
the second class constraints is 


X(S)ab^'= 

x{e)a 

X(w)ab* 

X{P)a^ 

18 

12 

18 

12 


which gives the total S = 60. Regarding the first class constraints, the situation is a 
little more complicated, due to the presence of Bianchi identities (see the Appendix). 
In particular, not all components of and are independent. To see this, 

take the derivative of to obtain 

. (45) 

The first term on the right-hand side is zero off-shell as a consequence of the second 
Bianchi identity (IA.8I) . The second term on the right-hand side is also zero off-shell, 
since it is a product of two constraints, 

i?^[“y 7r(R)/l*^' = 5'(i?)^[“y P(R)/l*^' = 0 . (46) 

Therefore, we have the off-shell identity 

V,(/)(i?)“^* = 0 , (47) 

which means that 6 components of (/)(i?)“^* are not independent of the others. In a 
similar fashion, we can calculate the following linear combination: 

V,(^(r)“ - = 

= - R<^\,ebk] (48) 

- i5(i?)-,,x(B)c*^' + . 

The term in the square brackets is zero off-shell as a consequence of the first Bianchi 
identity (IA.7I) . Additionally, the remaining two terms are products of constraints, and 
therefore also zero off-shell. Thus we have another off-shell identity, 

V.^(T)“ - \eM<l){Rr^^ = 0 , (49) 

which means that 4 components of (p{T)°-'' are not independent of the others. 

Taking (iTfl) and (|4^ into account, the number of independent components of the 
Hrst class constraints is 




ma 

4>{‘^)ab 

^!)(i?)“'’* 

Hor 

^(T)“ 

((.(VP)“'’ 

18 

4 

12 

6 

18-6 

4 

12-4 

6 


which gives the total of F = 70. Finally, substituting N, F and S into (Hll) . we obtain: 

n = 100 - 70 - y = 0 . (50) 

We conclude that the theory has no physical degrees of freedom. 
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As the final point of the analysis, we note that one can introduce the following 
canonical transformation on the phase space of the theory: 

7r(/3),*^ ^ + (1 - 20 e°^^^eak , 

/I \ n . 

7r(e)a* ^ 7r(e)a* = 7r(e)a* + ( 2 ~ W ’ 

while all other fields and momenta map identically onto themselves. It is easy to check 
that this change of variables is indeed canonical, since it does not change the Poisson 
structure. Moreover, the Hamiltonian (1361) and the primary and secondary constraints 
(|^ and (l40l) all transform such that 

(52) 

This is a symmetry of the action (11211 up to the boundary term, since 

S[1-0=S[0- [ e“A/3a. (53) 

J dM 

At the level of the full phase space, the canonical transformation dsa therefore maps 
between ^ and 1 — in particular between the BFCG theory (^ = 1) and the TPGT 
theory = 0). Nevertheless, after the elimination of the second class constraints 
and the phase space reduction, the situation will be more complicated, as we shall see 
in section O The canonical transformation between the BFCG and TPGT will still 
exist, but it will be singular in a certain sense, and not expressible in the generic form 
((?T|) . This will be discussed in detail in section [51 


4. Dirac brackets 


After the Hamiltonian analysis has been completed, we proceed to eliminate the second 
class constraints from the theory. This is done by introducing the Dirac brackets, 
defined as: 


{ F{t, x ), G{t, x') }d = { F(t, x ), G{t, x') } - 

f r (54) 

- d^y d^y'{F{t,x), x^it,y)}A/g{t,y,y'){x^{t,y'), G{t,x')} , 

JY. Jy 

where F and G are some functions of the phase space variables, while the kernel 
^( 43(^1 y, y') is the inverse of 

l:^^^{t,y,y’) = {x^(t,y), x^{t,y')} ■ (55) 


The multi-indices A and B count all 60 independent second class constraints. 

In order to evaluate the kernel and make the general definition (1541) more 
manageable, we proceed in several steps. First, from the Poisson brackets (l42l) we see 
that A^^{t,x,y) is diagonal in the space variables x and y, i.e. it can be written as 

A^^it,x,y) = A^^it)S^^\x-y). (56) 


That means that its inverse will also be diagonal in those variables, 
A;^B(t,y,y') = A^]j{t)6^^\y - y'), 

so that 


^d^y A^^{t,x,y)Aj^ly{t,y,y') = ScS^^\x - y '), 


(57) 

(58) 


Hamiltonian analysis of the BFCG theory for the Poincare 2-group 


13 


provided that 

= (59) 

From now on we will drop the explicit dependence of time from the notation of these 
matrices, for convenience. Substituting (l57l) into (l54l) and integrating over y', the 
Dirac brackets can be written in a simpler form 

{F,G}z5 = {F,G}-^ d^y{ F, x^{y) }A^^{ x^{y),G}, (60) 

where we have again simplified the notation by implicitly assuming appropriate 
spacetime dependence of variables F and G. 

As a second step, if we rewrite x^ and x^ as quadruples 

= (x(i3)“''*^x(w)“''\x(e)“,x(/3)“0 . 

yS = (x(i3)""’””,x(^)“"™,x(e)=’”,x(/3)=’"”) , 

we can write the matrix A"^^ in the block-diagonal form. 


(61) 


A^^ = 


/ ^abij\cdm 


\ 

^abi\cdmn 


V 

/\^aij\cm 

^ai\cmn 

/ 


(62) 


where we have used vertical bars to separate row from column indices, and the blank 
entries in the matrix are assumed to be zero by convention. According to (|42[) we have 


^abij\cdm _ ^a[c^d]h ^ 

^abi\cdmn _ _ ^^Oimn ^a[c ^d]b 

^ai\c 
/\aij\c 


■ \cmn 
'I cm 


_ _ Omni ac 

^ V 

_ ^Oijm^ac 


(63) 


The inverse matrix A^]^ then has a similar form, 

'' A"^ 

abij\cdm 


^AB — 


^ abi\cdmn 


V 


ai\cmn 


(64) 


aij\c 


Using this, from ([5^ one can obtain the equations 
1 


_ abij I cdm ~ 1 
2 ‘~^cdm\a'b' i' j' 

^abi I cdmn — 1 
^ ‘~^cdmn\a'b'i' 

az I cmn ~ 1 
2 cmn\a'i' 

A“^'l^™A"i, 

cm\a I j 




1 


(65) 


= 2<5“ , 

and then using (1631) one can solve them to obtain the components of the inverse matrix. 


A-i 


A-i 


by |cdm “ £OijmVa[cVd]b j 

abi\cdmn ~ ~^0imnValcVd]b j 

A“| 

azjcm? 


— —^0 


imnVac 5 


(66) 


A 


aij\cm 


— ^dijm'^O'C ■ 
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Here we have defined £0123 = — 

Finally, the third step is to substitute the matrix into (ISOl) in order to obtain 
an explicit expression for the Dirac brackets: 

{F,G}d = {F,G}- i£o.,fc [ d^y K^^'^iF, G, y) , (67) 

J 2 

where the kernel Kd^[F^ G, y) is 

K^^HF,G,y) = ^{F,xiBr^^^iy)}{xico)ab\y),G} 
-\{F,xiu^r\y)}{xiB)ab^’^iy),G} 

( 68 ) 

-{F,x{er{y)}{x{l3)a^'^{y),G} 

+{F,xiPr^iy)}{x{e)Jiy),G}. 

Having constructed the Dirac brackets, the next task is to express the constraint 
algebra in terms of them. This has two main consequences. The first is that the Dirac 
bracket between any quantity and any second class constraint is automatically zero, 
by construction. This is obvious from the definition (IMl) . The second is that after 
passing from Poisson brackets to Dirac brackets, the second class constraints can be 
set equal to zero off-shell, giving rise to the reduction of the phase space to one of its 
hypersurfaces. We will now concentrate on the constraint algebra, while the reduction 
of the phase space will be discussed in detail in the next section. 

Looking at the algebra of constraints, (HTj) . and (H51) . we see that it has the 
following rough structure: 

{(j), (j)} = (j), {x,x} = A, {(j)^x} = X, 

= {x,x} = 0, {0,x} = O. 

Knowing that the Dirac brackets between an arbitrary quantity and a second class 
constraint is zero by construction, we can immediately conclude that 

{x,x}D = d, {<?i',x}r> = 0, (70) 

which leaves only {^, (j)}D to be discussed. For this, a schematic calculation gives: 

{(?!),(/)}d = j {( t ), x }^~^{ x , 4 )} 

= + j xA"^x (71) 

= {0, 


due to the fact that the product of two constraints is zero off-shell. This is actually 
a proof that the algebra of primary constraints does not change when we pass from 
Poisson brackets to Dirac brackets. Therefore, we have: 


{cf{TY\ (fiVBUfo 

4>{VB),d}D 

{cf{VBY\ (f{VB),d}D 






(72) 
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while all other { 0, 4>}d vanish. 

As a final point, note also that for an arbitrary quantity A we have: 

{A,Ht}d = {A,HT}-j{A,x}^~^{x^HT} 

= {A,Ht}- J{A,x}^-^{x,^} (73) 

= {A,Ht}-J{a,x}a-\, 

where we have used the fact that the total Hamiltonian (1361) is a linear combination 
of first class constraints. The result can be rewritten as 

A = {A,Ht}d + J{A,x}^-^X, (74) 

which becomes the standard-looking equation of motion 

A = {A,Ht}d (75) 

when one reduces the phase space by promoting y 0 to off-shell equalities y = 0. 

This reduction is the subject of the next section. 

5. Phase space reduction 


The purpose of introducing Dirac brackets is to remove the second class constraints 
from the theory. When we use exclusively Dirac brackets, no result depends on second 
class constraints, and we can project all phase space points to the hypersurface defined 
by strong equalities y = 0, reducing its dimension from 2N to 2N — S, without 
changing any physical property of the theory, in particular without breaking its gauge 
symmetry. In our case, from (I40p we have the following off-shell equations: 

7r(H),bJ'= 

Tr(uj)ab" - £°"^^Babjk 

7r(/3)„y 

We will discuss these equations in two steps, first by analyzing the two equations 
independent of and then discussing the ^-dependent equations, which are more 
complicated. The first two equations can be rewritten as: 

n{B)aA^ = 0, B-\, = . (77) 

Note that we have expressed two conjugate variables in terms of other variables in 
the phase space. This reduces the full phase space to a hypersurface defined by 
these equations, namely orthogonal to the directions of TT{B)aG^ and diagonal in 
the directions of corresponding (i?,7r(a;)) planes. Given that we have eliminated 36 
phase space variables, the dimension of the hypersurface is 200 — 36 = 164, since 
the dimension of the full phase space was 2N = 200. On this hypersurface the 


0 , 

0 , 

0 , 

0 . 


(76) 
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expressions for some of the first class constraints (I39p simplify. In particular, the 
first four constraints remain unaffected, while the final four constraints become: 


cfiGY = + VMer , 

= Vi7r(a;)“''* - 2e[“i7r(e)''l* - . 


(78) 


Let us now turn to ^-dependent equations in (l76l) . We immediately see that there 
are two mutually incompatible cases, namely ^ 7 ^ 0 and C 7 ^ 1 - In the ^ 7 ^ 0 case, we 
see that we can solve the equations for e and 7r(e), 

e“i = , 7r(e)a* = , (79) 

again expressing two conjugate variables in terms of remaining ones. This reduces 
the dimension of the hypersurface even further, down to 164 — 24 = 140. If one 
does not impose any gauge fixing in the theory, this is the minimal dimension of the 
hypersurface, since we have in total 5" = 60 second class constraints. The final three 
first class constraints simplify even further, as follows: 

cf{GY = , 

HTr = -^V,7^(/3)“^ (80) 

cfiVBY^ = V.^(w)“''* - . 

We should stress that none of these equations make sense in the case ^ = 0, i.e. for the 
topological Poincare gauge theory (HU, but are completely valid for the case ^ = 1, 
which represents the BFGG theory ©• 

Alternatively, in the ^ 7 ^ 1 case, it is not a good idea to solve (|75)) for e, 7 r(e), 
since this cannot be done if ^ = 0. Instead, we can solve for /3 and 7 r(/ 3 ), 

= J^eo^Jk7T{eY '^, TT{f3)Y^ = ■ (81) 

This time the phase space reduces to another hypersurface, different from the previous 
one, but again of the same dimension 164 — 24 = 140. The final three first class 
constraints simplify again, however not to (1501) . but to: 


HGT 



(j){TY" 


(82) 

(f{VBY^ 




In this case the choice ^ = 1 does not make sense, which means that the BFGG 
theory case is excluded. Nevertheless, the case ^ = 0 is included, describing topological 
Poincare gauge theory. 

It is interesting to ask what happens in the case of generic when it is neither 
zero nor one. In that case one can solve (l76)) either for (e,7r(e)) or for (/3,7r(/3)). The 
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constraints can be expressed in either form (1501) or (I50|) . It is important to note, 
though, that the resulting hypersurface depends on the choice of In this case one 
can calculate the Dirac brackets 


{efy,7r(e),^}z5 = (I-aW^'^ 

(83) 


(84) 


It is then easy to verify that m is a canonical transformation from (e,7r(e)) to 
(/3, with (1701) being the inverse transformation. In particular, as long as ^ 7 ^ 0,1, 

this transformation maps (150|) to (ISO)) , and in addition maps (1551) to (1501) . justifying 
its canonical nature. 

However, the cases ^ = 0 and ^ = 1 are singular, and the canonical transformation 
dEH), 0701) does not make sense for either of those. Nevertheless, there exists a singular 
canonical transformation which maps between those two cases (see [H]), given as: 

= -eo^jMer ^, 7TiP)a^^ = . (85) 

In particular, for the case ^ = 0 the Dirac bracket (1551) evaluates to 

{e<^^,7rieW}D=S^SiS^^\ ( 86 ) 

while the constraints (1551) become 



= VMeY\ 




(87) 

(f{VBY^ 

= V,7r(a;)“^* - 2e[fy7r(e)''l*. 



On the other hand, when ^ = 1 we have from ((84l) 

{ , 7r(/3)r" }d = , (88) 

and from (l80l) we obtain 

mr = , 

HTr = -V,^(/3)“fy (89) 

= V,7r(a;)“''* - . 

The transformation (l85l) then maps (l86l) to (1551) and (1571) to (15^ . with its inverse 
mapping everything the other way around. 

To sum up, we have the following general situation. For a generic ^ one can write 
the theory using either (e, 7r(e)) variables or (/3,7r(/3)) variables, and there is a canonical 
transformation connecting these two sets of variables, for the same value of fy For 
the singular cases ^ = 0 and ^ = I one does not have a choice which variables to use, 
but there exists the canonical transformation (l85)) which maps the ^ = 0 theory into the 
^ = I theory, and vice versa. This canonical transformation is called singular because 
it cannot be obtained as a solution of the second class constraints (1751) . In contrast 
to ii]), which maps between various variables on the same hypersurface determined 
by the choice of fy the singular canonical transformation maps between two different 
hypersurfaces determined by choices ^ = 0 and ^ = I. This establishes the relationship 
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between the canonical structure of the BFCG model and the canonical structure of 
the topological Poincare gauge theory. 

We should also discuss the status of the canonical transformation (EU, which 
maps the full phase space onto itself, while inducing the transformation ^ —>■ 1 — 
Note that (ISTj) assumes that both momenta 7 r(e)a* and 7 r(/l)a*-’ are variables in the 
phase space. However, after the reduction of the phase space using either (IT^ or 
m, one of the two momenta is eliminated, and is not a variable in the reduced phase 
space. Therefore, one cannot formulate the canonical transformation (151|) as it stands, 
on the reduced phase space. Geometrically, every choice of ^ specifies one particular 
reduced phase space, and the symmetry ^ ^ 1 — ^ now maps between different spaces. 
While it is possible to construct a set of canonical transformations analogous to dSlD, 
in the sense of the map C 1 — these are not derivable from (I5ip . In particular, in 
the case of reduction (17^ , the change of variables 




(90) 


implements the map £, —)■ 1 — ^ in the constraints (1501) . Similarly, in the case of 
reduction (I5T]) . the change of variables 



(91) 


implements the same map in the constraints (1821) . Note that neither (1001) nor (10T|) 
can be derived from (ED), and also that both transformations are defined only for 
^ 7 ^: 0,1. Finally, for the case ^ = 0 —>• ^ = 1 and vice versa, we have the canonical 
transformation (l85l) . which also cannot be inferred from ED- We can of course infer 
the existence of all these canonical transformations from the fact that the map ^ ^ 1 —^ 
is a symmetry of the theory, see (I53|) . But the actual forms of the transformations 
cannot be obtained from each other, due to the different sets of variables in respective 
phase spaces. 

6. Conclusions 

Let us summarize the results. We have studied the Hamiltonian structure of the 
BFCG model ([D, and its relationship to the topological Poincare gauge theory (1TT1) . 
In section!^ we have defined both theories, proved the equivalence of their respective 
actions and Lagrange equations of motion, and introduced a generalized action ED 
which depends on a real parameter ^ and which reduces to the BFCG theory for ^ = 1 , 
while it reduces to the topological Poincare gauge theory for ^ = 0. In this way, we 
could perform the Hamiltonian analysis of both theories simultaneously. This was done 
in section [3l where it was established that there are no physical degrees of freedom, 
equation ED, the algebra of constraints (ED, ED and (l43|) has been computed, 
and the Hamiltonian of the theory written as a linear combination of first class 
constraints, equation (l36)) . In section |4] we have constructed the Dirac brackets (l67l) . 
which facilitate the elimination of the second class constraints from the theory without 
breaking its gauge symmetry. The geometrical consequences of this were explored in 
section [D The elimination of second class constraints projects the phase space onto 
one of its hypersurfaces, depending on the choice of the parameter Only at this 
point the difference between the BFCG theory and the topological Poincare gauge 
theory becomes observable, since they live on distinct hypersurfaces. For the generic 
hypersurface, when ^ 7 ^ 0 , 1 , the second class constraints can be solved in terms of 
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different sets of variables, and for every choice of ^ there is a canonical transformation 
which maps between these, mapping the hypersurface ^ onto itself. However, this was 
not possible for the singular cases ^ = 0,1. Instead, in these two cases there exists 
a singular canonical transformation which maps the ^ = 0 hypersurface to ^ = 1 
hypersurface and vice versa, establishing the equivalence between the BFCG theory 
and the topological Poincare gauge theory, and clarifying the relationship between 
their respective variables. 

The results obtained in this paper represent the straightforward generalization of 
results obtained in m of the BF theory based on the Lorentz group to the BFCG 
theory based on the Poincare 2-group. Some of the material presented overlaps with 
[12], but also improves on those previous results, in the following important ways. 
First, in this paper we have performed the full Hamiltonian analysis, as opposed to 
the shorthand procedure used in m- This facilitates a better basis for the future 
analysis of the constrained BFCG theory (|3|), whose relevance is very high since it 
is equivalent to GR. Second, the procedure used in [T2] was performed by employing 
a partial gauge fixing. This makes the calculations much simpler, but prevents us 
from computing the full algebra of constraints, in particular (1421) and (l4^ . In this 
paper the gauge symmetry was kept intact, the full algebra of constraints has been 
computed, and proved to be closed. And third, in this paper we have given a more 
detailed analysis of the relationship between the BFCG model and the topological 
Poincare gauge theory, providing a better insight into the geometry of the reduced 
phase spaces for both theories. 

The quantization of the theory has not been discussed, for two reasons. First, 
it was previously discussed in [12]. Second, while the formalism of Dirac brackets 
lends itself nicely to the canonical quantization, the quantum version of the BFCG 
model itself is of limited interest. Instead, the relevant theory is the constrained 
BFCG model ([3]), whose canonical quantization should give rise to a realistic theory 
of quantum gravity. This paper represents a first step toward the construction of 
such a theory, but the quantization procedure is however out of its scope, and will be 
studied in future work. 

Another topic for future work would be the Hamiltonian analysis of the BFCG 
model for an arbitrary 2-group [ido]. On one hand, models based on the extensions 
of the Poincare 2-group can give rise to a theory of quantum gravity with matter fields, 
providing a basis for the study of unification of gravity with other interactions at the 
quantum level. On the other hand, various choices of 2-groups can prove fruitful in 
the construction of new topological invariants of manifolds, which is a very important 
area of investigation in modern mathematics. For example, a state sum model based 
on the Euclidean 2-group has been constructed in [14] (see also [15]). One could 
therefore study the Hamiltonian structure of that model and perform the canonical 
quantization. Its comparison to the state sum model should provide greater insight 
into the structure of the theory. Various other applications may also be possible. 
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Appendix. Bianchi identities 


Recalling the definitions of the torsion and curvature 2-forms, 

= de“-h A e^ 7?“''= du;“^-h A (A.l) 

one can take the exterior derivative of T°- and i?“, and use the property dd = 0 to 
obtain the following two identities: 

Vr“ = dT“-h a;“b A T'’ = A e^ .. . 

= di?“'’ -h a;“c A -h A = 0 . 

These two identities are universally valid for torsion and curvature, and are called 
Bianchi identities. By expanding all quantities into components as 

= ^TV^da:'^ A dx" , A da:" , (A.3) 


e“ = e^dx'^, , (A.4) 

and using the formula dx^ A dx" A dx'’ A dx'^ = e^^'^'^d^x, one can rewrite the Bianchi 
identities in component form as 

eV-P - R\^,e\) = 0, (A.5) 

and 


= 0. (A.6) 

For the purpose of Hamiltonian analysis, one can split the Bianchi identities into 
those which do not feature a time derivative and those that do. The time-independent 
pieces are obtained by taking A = 0 components: 

- R\,,e\) = 0 , (A.7) 


^Oijky^Rab^^ = 0 . 


(A.8) 


These identities are valid as off-shell, strong equalities for every spacelike slice in 
spacetime, and can be enforced in all calculations involving the Hamiltonian analysis. 
The time-dependent pieces are obtained by taking X = i components: 

^o^3k (Vor“,fe - 2V,T%fc - 2R\o,e\ - R\jke\) = 0 , (A.9) 

and 

(VoR“^fe - 2V,R“%fc) = 0. (A.IO) 

Due to the fact that they connect geometries of different spacelike slices in spacetime, 
they cannot be enforced off-shell. Instead, they can be derived from the Hamiltonian 
equations of motion of the theory. 

In light of the Bianchi identities, we should note that the action ([7]) features two 
more fields, /3“ and which also have field strengths G“ and VH“^, and for which 
one can similarly derive Bianchi-like identities, 

VG“ = R“^,A/3^ = R“c A 5'='’-h A (A.ll) 
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However, due to the fact that both fi°‘ and are two-forms, in 4-dimensional 
spacetime these identities will be single-component equations, with no free spacetime 
indices, 

)=0, (A.12) 

and similarly for V^B°‘^. Therefore, these equations necessarily feature time 
derivatives of the fields, and do not have a purely spatial counterpart to (IA.7I) and 
(|A.8I) . In this sense, like the time-dependent pieces of the Bianchi identities, they do 
not enforce any restrictions in the sense of the Hamiltonian analysis, but can instead 
be derived from the equations of motion and expressions for the Lagrange multipliers. 
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